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Q\ ' Abstract 

o 1 + 

In this article, we calculate the masses and the pole residues of the ^ 
heavy baryons Q®(css) and Q^(bss) with the QCD sum rules. The numerical 
X> ■ values M n o = (2.72 ± 0.18) GeV (or M n o = (2.71 ± 0.18) GeV) and M n - = 

rj^ ' (6.13 ± 0.12) GeV (or M Q - = (6.18 ± 0.13) GeV) are in good agreement with 

^ ■ the experimental data. 

PACS number: 14.20.Lq, 14.20.Mr 
r-| . Key words: Q®(css), Q^(bss), QCD sum rules 
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1 Introduction 

J, . 

The charmed and bottomed baryons which contain a heavy quark and two light 
j> . quarks are particularly interesting for studying dynamics of the light quarks in the 

presence of a heavy quark. They behave as the QCD analogue of the familiar 
hydrogen bounded by the electromagnetic interaction, and serve as an excellent 
ground for testing predictions of the constituent quark models and heavy quark 
symmetry (TJ [2]. The u, d and s quarks form an SU(3) flavor triplet, 3x3 = 3 + 6, 
two light quarks can form diquarks of a symmetric sextet and an antisymmetric 
O ! antitriplet. For the S'-wave baryons, the sextet contains both spin-| and spin-| 

j> \ states, while the antitriplet contains only spin-| states. 

There have been many theoretical approaches dealing with the heavy baryons 
with one heavy quark, for example, the relativized potential quark model [3], the 
Feynman-Hellmann theorem and semi-empirical mass formulas [I], the Faddeev 
equation [51 E], the combined expansion in 1/mg and 1/N C [7], the lattice simu- 
lations [SI H], the relativistic quark-diquark approximation [TUl fTT] , the hyperfme 
interactions [12J, the QCD sum rules in the heavy quark effective field theory 
the variational approach [H], and the full QCD sum rules [151 H], etc. 
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The | + antitriplet states (A+, H+, H°), and the ^ + and | + sextet states (fl c , S c , E' c ) 
and (f2*,S*,S^*) have been established; while the corresponding bottomed baryons 
are far from complete, only the A b , E b , ££, E b have been observed [T7j. 

Recently the DO collaboration reported the observation of the doubly strange 
baryon in the decay channel Q b — > J/ipQ~ (with J/ip —>■ fi + n~ and Q~ — > 
AK~ — > pn~K~) in pp collisions at ^fs = 1.96 TeV [IB]. The experimental value 
M n - = 6.165 ± 0.010 ± 0.013 GeV is about 0.1 GeV larger than the existing theo- 
retical calculations [H El [71 El El [101 HH [12l UHl [H], while the theoretical prediction 



1 E-mail, wangzgyiti@yahoo.com.cn. 



1 



M n o « 2.7 GeV [H El d El [91 HQl HH H21 [131 HH is consistent with the experimental 
data M n o = (2.6975 ± 0.0026) GeV [Tf] . 

In previous work, we have calculated the masses and the pole residues of the | 
heavy baryons f2* and Ql with the QCD sum rules [12] . In this article, we extend 
our previous work to study the corresponding ^ + heavy baryons and Q b . 

In Ref . |20] . the masses of the heavy baryons and are studied in the full 
QCD sum rules with the same interpolating current as the present work. The authors 
choose the tensor structure p 1 and obtain different spectral density 0. 

The masses of the Aq, Eq and 5q have been calculated with the full QCD sum 
rules [H3, HEl [20]. The masses of the Eq, Eq and Aq have been calculated with the 
QCD sum rules in the leading order of the heavy quark effective theory [2U [221 [23], 
and later the 1/rriQ corrections were studied [2U (25J [26]. Furthermore, the masses 
of the orbitally excited heavy baryons with the leading order approximation [271 12H] 
and the l/rnq corrections [29J in the heavy quark effective theory have also been 
analyzed. Recently the | and | bottomed baryons were studied with the QCD 
sum rules in the heavy quark effective theory including the l/rnq corrections |13j . 

In the QCD sum rules, the operator product expansion is used to expand the 
time-ordered currents into a series of quark and gluon condensates which parame- 
terize the long distance properties of the QCD vacuum. Based on the quark-hadron 
duality, we can obtain copious information about the hadronic parameters at the 
phenomenological side [30], EI] . 

The article is arranged as follows: we derive the QCD sum rules for the masses 
and the pole residues of the and fl^ in section 2; in section 3 numerical results 
are given and discussed, and section 4 is reserved for conclusion. 



2 QCD sum rules for the Q c and Q b 

In the following, we write down the two-point correlation functions U a (p) in the 
QCD sum rules approach, 

J*( X ) = e^ k S J(x)C^ Sj (x) l5 rQt(x), (2) 
X a N(p,s) = (0\J a (0)\Q a (p,s)), (3) 

where the upper index a represents the c and b quarks respectively, i, j and k are 
color indexes, C is the charge conjunction matrix, the N(p, s) and A a stand for the 
Dirac spin vector and the pole residue of the heavy baryon Q a , respectively. 

2 Whcn I finish the work and submit it to the network |http: / /arXiv.org] I learn from Prof. 
Nielsen that they have studied the masses of the heavy baryons 51" and £1^ in the full QCD sum 
rules with the same interpolating current J a (x) [20] . 
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The correlation functions U a (p) can be decomposed as follows: 

n<» = *n»(p) + n»(p) , (4) 

due to Lorentz covariance. The first structure ^ has an odd number of 7-matrices 
and conserves chirality, the second structure 1 has an even number of 7-matrices and 
violate chirality. In the original QCD sum rules analysis of the nucleon masses and 
magnetic moments [321 E21 EH E2] , the interval of dimensions (of the condensates) 
for the odd structure is larger than the interval of dimensions for the even structure, 
one may expect a better accuracy of the results obtained from the sum rules with 
the odd structure. 

In this article, we choose the two tensor structures to study the masses and the 
pole residues of the heavy baryons Q° and f^~, as the masses of the heavy quarks 
break the chiral symmetry explicitly. 

The components II* (p) with the odd structure have smaller dimension of mass 
than the components II^Qo) with the even structure, naively we expect the compo- 
nents IIJ(p) have better convergent behavior in the operator product expansion, and 
the numerical results confirm this conjecture, see Tables 1-2. 

We insert a complete set of intermediate baryon states with the same quantum 
numbers as the current operators J a (x) into the correlation functions U a (p) to obtain 
the hadronic representation [301 EI]- After isolating the pole terms of the lowest 
states fl a , we obtain the following result: 

n°M = a^±* + .... (5) 

Ivl n a P 

In the following, we briefly outline the operator product expansion for the corre- 
lation functions U a (p) in perturbative QCD. The calculations are performed at large 
space-like momentum region p 2 <C 0, which corresponds to small distance x ~ re- 
quired by validity of the operator product expansion. We write down the "full" 
propagators Sij(x) and Sq(x) of a massive quark in the presence of the vacuum 
condensates firstly [31] . 



1152 y y 32vr 2 x 2 ^ 



S l A(x) = — — / d A ke 



a/3 

(2tt) 4 \j6-m Q 4 (k 2 - »•- - 



* I .1. -ft-J <% &Aj Vaf3$ + m Q ) + (ft + m Q )a a0 



m, 



7T 2 .a s GG. h 2 + m Q f{ . 

+ i ( ^ m °(^ + '^' (6) 

where (sg s aGs) = (sg s a af sG a ^s) and ( QsGQ ) = ( asGa ^ G — ^ then contract the quark 
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fields in the correlation functions U a (p) with Wick theorem, and obtain the result: 

U a (p) = -2ie ijk e i,j ' k ' J d i xe i ^Tr[ 1 ^x) lu CSl l {x)C] l5l ,S k Q k \x) lvl ,.{l) 

Substitute the full s, c and b quark propagators into above correlation functions and 
complete the integral in the coordinate space, then integrate over the variable k, 
we can obtain the correlation functions II" (p) at the level of quark-gluon degree of 
freedom. Once analytical results are obtained, then we can take the quark-hadron 
duality and perform Borel transform with respect to the variable P 2 = —p 2 , finally 
we obtain the following two sum rules with respect to the tensor structures p 1 and 1 
respectively: 

M 2 1 C s ° a f 1 

A a ex P(-^) = 1(^4/ ds J^dxx(l-xf {m 2 a - s) (3m 2 a - 5 s) exp( 

^ m s (ss) f ds f dxx(2 - 3x)exp(--^) 
71 Jth JA a M 
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1 ,a s GG, f s ° , f 1 , , A r , , s v 



IT 2 
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6tt 2 J y M 2 Ly M 2 

1 ,a s GG. f 1 , ,2a; 3 - 9x 2 + 9x - 2 (l-a;) 3 m 2 
(— ^" ) / d < Z ZTH ) 



1447T 2 7T Jo X XM 2 

~2 / ™n 

m a exp 



M 2 

2{ss) 2 , m 2 m s {sg s aGs) , m 2 



^-W^^J^^<-jf 2 ), (8) 



/ M l ^ 



«1 



3m a m, 



96tt 2V rr 



'~ ] L ds J Aa dxi ~ 2 - 2x - 3)eM -^ : 

a s GG f 1 (l-x) 2 „ 2 m 2 



96tt 2 ^ tt ; 7 x- 

5m a m s (sg s aGs) m\ Am a (ss) 2 m 2 

+ 12^ eXp( " + 6XP( - ' 

(9) 
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where th = (m a + 2m s ) 2 , A a = — and fh 2 = —. 

We carry out the operator product expansion to the vacuum condensates adding 
up to dimension-6 and calculate the same Feynman diagrams as in the sum rules for 
the | + heavy baryons [T9j (where the two tensor structures g^ v $ and are chosen), 
the contribution from each of the Feynman diagrams differs from the corresponding 
one in our previous work. 

In calculation, we take assumption of vacuum saturation for the high dimension 
vacuum condensates, they are always factorized to lower condensates with vacuum 
saturation in the QCD sum rules, factorization works well in large N c limit. In 
this article, we take into account the contributions from the quark condensates 
(ss), (ss) 2 , mixed condensate (sg s aGs), gluon condensate ( a " ), and neglect the 
contributions from other high dimension condensates, which are suppressed by large 
denominators and would not play significant roles. 

Differentiate the above sum rules with respect to the variable -rk-, then eliminate 
the pole residue A^ a , we obtain two QCD sum rules for the masses M^ a with respect 
to the tensor structures^ and 1, respectively. 

In the heavy quark limit tuq — > oo, the quarks c and b degenerate, we replace 
the heavy quark field Q a (x) in the baryon current J a (x) with the effective field 
h(x) to obtain the interpolating current Jh(x) = e l: > k sj (x)C^ ^s j(x)^^hk(x) in the 
heavy quark effective theory, which differs from the corresponding one constructed 
from the heavy quark effective theory remarkably and does not warrant the tensor 
structure ^ [22]. 

In the heavy quark limit, the correlation function IT a (u;) (with & — v ■ p) can be 
written as IP (a;) =^Ui(lj)+I12(uj)(I1i(u)) ^ U. 2 (uj)). In Ref . |13] . the current operator 
r)(x) = e^ k sj {x)C^^Sj{x)^^[^— /)v^}hk(x) is used to interpolate the baryons Q a , 
and the correlation function is written as U a (uj) = ^tl(uj). 

The effective heavy quark propagator has the following form, 

S a h b (x) = 5 ab ^ !°° dtS(x-vt), (10) 
^ Jo 

where the t> M is a four- vector with v 2 — 1. The calculations of the operator product 
expansion can be performed in the coordinate space and are greatly facilitated, as 
we do not need the mixed picture both in coordinate and momentum spaces. We 
calculate the same Feynman diagrams as in the full QCD, the contributions of some 
diagrams vanish in the heavy quark limit. 

Finally we obtain two sum rules with respect to the tensor structures ft and 1 
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respectively: 

^9 , A. 1 f^ , , , uj . m s (sq s aGs) (ss) 2 

F exp(— — ) = - ^exp(--) + + i_L , (11 ) 

/ A\ 1 P j 5 / 3m s (ss) r° . , u. 

F exp -- = -— / dujuj exp -- — / dujuexp -- 

T lOvr 4 J 2ms T vr 2 J 2ms T 

1 ,a s GG. r° , . cj x 5m s (sg s aGs) 2(ss) 2 

7 (— )/ rfcjwexpf ) + -+ x 7 , 

16vr 2V tt ' J 2ms l K T ! 24tt 2 3 ' 

(12) 

where we have used the definitions for the binding energy A and the pole residue F, 
UmitmQ^ooMn = rriQ + A and (0| J(0)\fl(v , s)) = FN(v,s), the N(v,s) is the Dirac 
spin vector in the heavy quark limit. 

Differentiate the above sum rules with respect to the variable ^ then eliminate 
the pole residue F, we obtain two QCD sum rules for the binding energy A with 
respect to the tensor structures jj and 1, respectively. 



3 Numerical results and discussions 

The input parameters are taken to be the standard values (qq) = — (0.24±0.01 GeV) 3 , 
(ss) = (0.8 ± 0.2) (qq), (sg s aGs) = m 2 (ss), m 2 = (0.8 ± 0.2) GeV 2 , (^p) = 
(0.33 GeV) 4 , m s = (0.14±0.01) GeV, m c = (1.4±0.1) GeV and m b = (4.8±0.1) GeV 
at the energy scale about \i = 1 GeV [30, 31, 36J. The contribution from the gluon 
condensate ( Qs ) is less than 3%, and the uncertainty is neglected here. 

We usually consult the experimental data in choosing the Borel parameter M 2 
and the threshold parameter sq. There lack experimental data for the phenomeno- 
logical hadronic spectral densities of the bottomed baryons at present, only the 
Eft, E£, E b [T7j and [TS] have been observed, we can borrow some ideas from the 
light baryon spectra [TT] . 

For the octet baryons with J p = | , the mass of the proton is M p = 938 MeV, 
and the mass of the first radial excited state A^(1440) (the Roper resonance) is 
M 1440 = (1420 - 1470) MeV w 1440 MeV; the mass of the ground state E is M s = 
1189.37 MeV, and the mass of the first radial excited state E(1660) is M 1660 = 
(1630-1690) MeV^ 1660 MeV; the mass of the ground state H is Ms = 1321. 7MeV, 
while the spin-parity of the high excited states 2(1620), 2(1690), 5(1950), 5(2030), 
5(2120), 5(2250), 5(2370), 5(2500) have not been determined yet pi]. For the 
decuplet baryons with J p = | + , the mass of the A(1232) is M1232 = (1231 — 
1233) MeV w 1232 MeV, and the mass of the first radial excited state A(1600) is 
Mieoo = (1550 - 1700) MeV w 1600 MeV; the mass of the ground state E(1385) 
is M1385 = 1382.8 MeV, and the mass of the first radial excited state E(1840) is 
Mi 840 = 1840 MeV; the ground states 5(1530), 0(1672) are well established, while 
the spin-parity of the high excited states £7(2250), fi(2380), Q(2470) have not been 
determined yet [TT] . 
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Eq.(8) 


Eq.(9) 


perturbative term 


+93% 


+70% 


(88) 


-1% 


+21% 


(sg s aGs) 


-1% 


-3% 


(ss) 2 


+8% 


+10% 


(^) 


+ 1% 


+ 1% 



Table 1: The contributions from different terms in the sum rules for the with 
the central values of the input parameters " Set I" . 





Eq.(8) 


Eq.(9) 


perturbative term 


+94% 


+76% 


(88) 


-2% 


+23% 


(sg s aGs) 


-1% 


-2% 


(ss) 2 


+8% 


+5% 


<^> 


+1% 


-3% 



Table 2: The contributions from different terms in the sum rules for the Q b with 
the central values of the input parameters " Set I" . 

From the experimental data for the baryons consist of qqq and qqs, we can see 
that the energy gap between the ground states and the first radial excited states is 
about (0.4 — 0.5) GeV, the SU(3) breaking effects are rather small. So in the QCD 
sum rules for the baryons with the light quarks, the threshold parameters sq are 
always chosen to be y/so = M gr + 0.5 GeV [32], [33] Ell EH], here gr stands for the 
ground states. The central values of the threshold parameters for the heavy baryons 
tt° c and tt; can be chosen as s° Qc = (2.6975+0.5) 2 GeV 2 and s^ b = (6.165+0.5) 2 GeV 2 , 
respectively. 

In this article, the threshold parameters and the Borel parameters are taken as 
s o c = (10.5 ± 1.0) GeV 2 and M 2 = (2.2 - 3.2) GeV 2 for the charmed baryon and 
s o^ = (44.5 ± 1.0) GeV 2 and M 2 = (5.0 - 6.0) GeV 2 for the bottomed baryon tt b ; 
thereafter these parameters will be denoted as "Set I". In the heavy quark limit 
rriQ — > oo, the corresponding Borel parameter T and threshold parameter ujq can be 
taken as T = (0.4 — 0.6) GeV and ujq = (1.7 — 1.9) GeV; and we will refer them as 
"Set III". 

We can make another estimation for the threshold parameters by adding the 
heavy quark masses to the ground state | + baryons which consist of qss, Sq c = 
(Ms + m c + 0.5 GeV) 2 « 10.5 GeV 2 and s° Q}> = (M s + m b + 0.5 GeV) 2 w 44.0 GeV 2 , 
which are almost the same values as " Set I" . 

The contributions from different terms at the central values of the input parame- 
ters are presented in Table. 1 and Table. 2, respectively. From the two tables, we can 
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n 6 - 


i> k, Set I 


(34 - 60)% 


(35 - 49)% 


1 & Set I 


(34 - 62)% 


(36 - 51)% 


i> & Set II 


(52 - 75)% 


(51 - 66)% 


1 & Set II 


(53 - 78)% 


(51 - 67)% 



Table 3: The contributions of the pole terms from different sum rules, here we take 
the central values of the input parameters except for the Borel parameter M 2 . 

expect convergence of the operator product expansion. The components II" (p) with 
the odd structure have smaller dimension of mass than the components IT^p) with 
the even structure, see Eq.(4), it is not unexpected that the components II" (p) have 
better convergent behavior in the operator product expansion and result in better 
QCD sum rules, see Table 4. 

In the heavy quark limit, the contribution from the quark condensate (ss) van- 
ishes for the odd structure tf, while the condensate (ss) has a numerical coefficient 
^2 for the even structure 1, see Eqs. (11-12), which can explain the hierarchy appears 
in the operator product expansion naturally, see Tables 1-2. 

The contributions of the pole terms for different sum rules are shown in Table 3, 
from the table, we can see that the pole dominance condition is marginally satisfied 
for the parameters " Set I" . 

We can choose smaller Borel parameters M 2 or larger threshold parameters s° 
to enhance the contributions from the ground states. However, if we take larger 
threshold parameter s°, the contribution from the first radial excited state maybe 
included in; on the other hand, for smaller Borel parameter M 2 , the sum rules are 
not stable enough, the uncertainty with variation of the Borel parameter is large. 

In this article, we also present the results with smaller Borel parameters M 2 and 
larger threshold parameters s° a , s° Uc = (11.0 ± 1.0) GeV 2 and M 2 = (1.9 - 2.6) GeV 2 
for the charmed baryon and s& = (46.5 ± 1.0) GeV 2 and M 2 = (4.5-5.4) GeV 2 
for the bottomed baryon . These parameters (denoted as "Set II") satisfy both 
pole dominance (see Table 3) and convergence of the operator product expansion. 

Taking into account all uncertainties of the input parameters, finally we obtain 
the values of the masses and the pole residues of the heavy baryons 0° and 
which are shown in Figs. 1-4 (and Table 4), respectively. 

In the heavy quark limit, we obtain the binding energy A (see Fig. 5) 

A = 1.30±0.17GeV, 

A = 1.10 ±0.26 GeV, (13) 

for the tensor structures f> and 1 respectively. The corresponding masses are pre- 
sented in Table 4. 

The values M n o = (2.72 ± 0.18) GeV and M n o = (2.71 ± 0.18) GeV from the 
sum rules with the tensor structure p 1 are in good agreement with the experimental 
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Upper bound; 
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Upper bound; 
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Figure 1: The mass Mqo from the sum rules with different tensor structures and 
input parameters, A, B, C and D correspond to ^ & Set I, 1 & Set I, & Set II and 
1 & Set II, respectively. 





M n o(GeV) 


M n -(GeV) 


Ano(GeV^) 


V-(GeV a ) 


i> k, Set I 


2.72±0.f8 


6.13 ±0.12 


0.088 ±0.026 


0.118 ±0.026 


f & Set I 


2.66 ±0.19 


6.09 ±0.11 


0.064 ±0.017 


0.088 ±0.017 


i> & Set II 


2.71 ±0.18 


6.18 ±0.13 


0.090 ±0.027 


0.133 ±0.030 


1 & Set II 


2.63 ±0.20 


6.12 ±0.13 


0.065 ±0.017 


0.099 ±0.019 


i & Set III 


2.70 ±0.27 


6.10 ±0.27 






1 & Set III 


2.50 ±0.36 


5.90 ±0.36 






Exp data 


2.6975 ± 0.0026 


6.165 ±0.013 







Table 4: The masses and pole residues from the sum rules with different tensor 
structures and input parameters. 
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Figure 2: The mass M n - from the sum rules with different tensor structures and 
input parameters, A, B, C and D correspond to ^ & Set I, 1 & Set I, $ & Set II and 
1 & Set II, respectively. 
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Figure 3: The pole residue A^ c from the sum rules with different tensor structures 
and input parameters, A, B, C and D correspond to $ & Set I, 1 & Set I, $ k, Set II 
and 1 & Set II, respectively. 
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Figure 4: The pole residue Xn b from the sum rules with different tensor structures 
and input parameters, A, B, C and D correspond to $ & Set I, 1 & Set I, $ h Set II 
and 1 & Set II, respectively. 











Central value; 




: a 


Upper bound; 






Lower bound. 









0.40 



0.45 



0.50 

T [GeV] 



0.55 



0.60 



1.6 
1.4 
1.2 
> 1.0 

Oo.8 

^ 0.6 
0.4 
0.2 
0.0 













Central value; 




B 


Upper bound; 






- Lower bound. 











0.40 



0.45 



0.50 

T [GeV] 



0.55 



0.60 



Figure 5: The binding energy A with variation of the Borel parameter T, A and B 
correspond to the tensor structures f and 1, respectively. 



12 



data Mqo = (2.6975 ± 0.0026) GeV [T7], other theoretical predictions also indicate 
the value' is about 2.7 GeV [11 El El El El HOI HH H21 H31 HI]- The experimental 
value M n - = 6.165 ± 0.010 ± 0.013 GeV is about 0.1 GeV larger than the existing 
theoretical calculations [H El El El El HTJ1 HU H21 H31 [U| (including the QCD sum 
rules in the heavy quark effective theory with l/rriq corrections [13]), our predictions 
M n - = (6.13 ± 0.12) GeV and M a - = (6.18 ± 0.13) GeV based on the sum rules 
with the tensor structure^ are excellent. 

The parameters "Set I" satisfy the pole dominance criterion marginally, the 
parameters "Set II" have larger threshold parameters than that of the "Set I", 
which may take into account some contributions from the high resonances. More 
experimental data are needed to select the ideal sum rules. Once reasonable values 
of the pole residues Xn c and Xn b are obtained, we can take them as basic input 
parameters and study the hadronic processes with the (light-cone) QCD sum rules. 

In Ref. J2U] , the tensor structure jf is chosen, and the spectral density (including 
the perturbative term, quark condensate term (ss), mixed condensate term (sg s aGs) 
and gluon condensate term ( ° 3 )) differs from mine, the most evident difference 
is that the terms of the vacuum condensate (ss) present in Eq.(8) are absent. I 
check all calculations carefully and confirm my results. For example, although the 
term (ss) vanishes in the heavy quark limit (see Eq.(ll)), there indeed exist such 
a term in the full QCD. The threshold parameter s^ c = (10.0 — 11.5) GeV 2 and 
the prediction M^o = (2.65 ± 0.25) GeV are consistent with mine, while the lower 
bound of the threshold parameter = (41.0 — 45.0) GeV 2 is much smaller than 
mine, and the prediction M n - = (5.82 ± 0.23) GeV differs from the experimental 
data remarkably. 

The sum rules in the heavy quark limit with the odd structure f also result in 
excellent values M Q o = (2.70 ± 0.27) GeV and M a - = (6.10 ± 0.27) GeV, although 

c b 

the uncertainties are large; the values from the even structure 1 are not good enough, 
see Table 4, the 1/toq corrections maybe large. In Ref. [13], the threshold parameter 
is taken as u = 1.55 GeV, which differs from mine remarkably. The predictions 
(the central values are M n o = 2.62 GeV and M n - = 5.97 GeV without the I/thq 
corrections) are smaller than mine about (0.08 — 0.10) GeV. It is not unexpected, 
different interpolating currents and input parameters can lead to different results. 

In this article, we do not take into account the next-to-leading order corrections 
to the perturbative term, the corrections maybe large. In the massless limit m s = 
m u = rrid = 0, we can resort to analytical expressions of the correlation functions for 
the heavy baryon currents with one heavy quark in the finite mass limit [3T| I3"8"] or 
infinite mass limit [39, 40J to make possible improvement for the predicted masses 
Mq ci , while the next-to- leading order corrections to the perturbative terms for the 
baryon currents with three massless quarks are calculated in Refs. [HI H2J H3l HI] . 
The analytical expressions of the perturbative a s corrections presented in Ref. |37j 
are lengthy enough; while the expressions in the heavy quark limit mq — > oo are 
simpler and more easy to deal with, however, the interpolating current differs from 
mine remarkably [39~l HQ]. The mass of the s quark plays an important role, the 
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fi& and have equal masses in the limit m s = 0, from the experimental data 
M Ej = (5807.8 ± 2.7) MeV [TT] , we can see that the flavor SU(3) breaking effects 
(about 350 MeV) are rather large, m s = m u = = maybe a crude approximation. 
We prefer another work to study those perturbative effects in detail. 

4 Conclusion 

In this article, we employ the QCD sum rules to calculate the masses and the pole 
residues of the heavy baryons f2°(css) and Q^(bss) including the contributions of the 
vacuum condensates adding up to dimension six in the operator product expansion. 
The numerical values M n o = (2.72 ± 0.18) GeV (or Moo = (2.71 ± 0.18) GeV) and 
Mo- = (6.13 ±0.12) GeV (or M n - = (6.18 ±0.13) GeV) are in good agreement with 

b b 

the experimental data, while the existing theoretical predictions for the mass M n - 
is about 0.1 GeV lower than the experimental value. We can take the pole residues 
Aq c and Xn b as basic parameters and study the hadronic processes, for example, 
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